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Abstract 

Neutrino-nucleus elastic scattering is reduced in dense matter because of 
correlations between ions. The static structure factor for a plasma of electrons 
and ions is calculated from Monte Carlo simulations and parameterized with 
a least squares fit. Our results imply a large increase in the neutrino mean 
free path. This strongly limits the trapping of neutrinos in a supernova by 
coherent neutral current interactions. 
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A (core collapse) supernova radiates large numbers of neutrinos. Indeed, the energy 
in neutrinos is 100 times greater then that in all other forms of matter Jy]. Therefore, 
supernova models may depend on the details of neutrino interactions in dense matter. In 
this paper, we calculate how correlations in the medium modify the important neutrino- 
nucleus elastic scattering cross section. This cross section is large because it involves coherent 
scattering from all of the nucleons in a nucleus || . However, when the neutrino wave length 
is comparable to the inter particle spacing there are also coherent contributions from different 
nuclei. These can screen the interaction and lead to a large reduction in the cross section. 

In the present supernova model, the core of a massive star runs out of nuclear fuel 
and collapses [||. This core is composed of a dense plasma of electrons and nuclei. As 
the density reaches 10 11 to 10 12 g/cm 3 the medium starts to become opaque to neutrinos. 
The neutrino opacity is thought to be dominated by neutrino-nucleus elastic scattering (as 
long as a significant fraction of the matter is in nuclei). This opacity insures that neutrino 
transport involves diffusion (rather then free streaming). The diffusion time can become long 
compared to the dynamical time scale, thus trapping neutrinos and their lepton number. 

The neutrino-nucleus elastic cross section in free space is [|J, 

G 2 C 2 E 2 (l + cosfl) 
dao/dll = — 2 , (1) 

with G the Fermi constant, E v the neutrino energy, 9 the scattering angle and the total 
weak charge C of a nucleus of charge Z and neutron number N is 

C = -2Zsm 2 Q w + (Z- AO/2. (2) 

(We use a Weinberg angle of sin 2 Q\y = .223.) In a dense plasma this cross section is 
modified by electron ]5]|| and ion || screening. Imagine a single impurity ion in a dense 
plasma. Extra electrons will be attracted to the impurity. Since these electrons have weak 
interactions they screen both the electro-magnetic and weak charge of the ion. However, 
the very dense relativistic electron gas is quite rigid because of the large Fermi momentum. 
This limits the effect of electron screening (see below). 

Other ions can also screen the impurity by creating a small hole in the ion distribution. 
At temperatures of order one MeV, the ions are essentially classical and their screening is 
not impeded by a large Fermi energy. Therefore, we will focus on ion screening in this paper. 
Some ion screening results have been presented in Ref. ||. Here we calculate screening for 
a broad range of densities and determine its impact on the neutrino mean free path. We 
also provide a parameterization of our results. This, or a further simplification, will allow 
the incorporation of screening in neutrino transport codes. 

Ion screening is included by multiplying Eq. (fl|) by the static structure factor S q of the 
ions 0, 

da/dQ = da /dfl S q . (3) 

Here q is the momentum transfer and da/dfl the effective cross section in the medium. We 
discuss S q below. 

The transport cross section is the angle integral of Eq. (^) with a factor of (1 — cos#), 
o l = J da/dVl (1 - cos) dSl = a* < S > . (4) 
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The free transport cross section is, a l Q = %G 2 C 2 E 2 jTT ^ and < S > is the angle average of S q , 

<S>=j[ dcos9 (1 + cos0)(l - cos9) S q(e) . (5) 

Here (1 + cos8) is from the angular dependence of the free cross section and q(6) 2 = 2E 2 (1 — 
cos#). Thus, ion screening can be incorporated into neutrino transport codes be multiplying 
the existing interactions by the factor < S > []. This depends on the density, temperature 
and neutrino energy. The transport mean free path A then follows, A = l/(ncr*), with n the 
number density of ions. 

The static structure factor S q is determined from a Monte Carlo simulation || of the 
radial distribution function g (r) ||. This gives the probability to find another ion a distance r 
from a given ion and is calculated by histograming the relative distances in the simulation || , 

S q = 1 + n f d 3 re- ixi r {g(r) - 1). (6) 



Equations (|5]||) yield a simple integral for < S >, 



<S>=1 + — drf(2E„r)(g(r)-l), (7) 



o 

with 

f{x) = 72(cosx + xsinx — l)/x 4 — 6(5cosx + xsinx + l)/x 2 . (8) 

The classical canonical partition function is simulated using iVj w 1000 ions in a box 
of volume V = Ni/n with periodic boundary conditions. The ions interact via screened 
Coulomb potentials, 

Z 2 p 2 

v ( r ) = — e -r/x e _ ( 9 ) 



Here A e = tc/ (ekp) describes the electron screening of the ion-ion interaction Note, this 



Yukawa approximation can be replaced by a more accurate description at high momentum 
transfers. However, we are primarily interested in momentum transfers q much less then the 
Fermi momentum q « kp. Therefore Eq. (^J) should be adequate for our purposes. 

The system is warmed up for about 200 Metropolis sweeps starting from either a simple 
cubic lattice or a uniform distribution. Statistics are then accumulated using 400 config- 
urations each of which is separated by 5 sweeps. This yields S q with a typical statistical 



accuracy of 2 — 3 x 10 . These results are close to S q for a pure one component plasma [JIT 

We parameterize our Monte Carlo results for < S > as a function of two dimension-less 
variables. It is a strong function of 

E = E u a (10) 



1 Note, Eqs. (|l],|]) have ignored axial current contributions to the cross section. These may be 
significant when < S > is very small. 
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Here, the ion sphere radius a measures the average distance between ions 0, 

a = [3/(47rn)] 1/3 . (11) 

Next, < S > is a weak function of V which characterizes the strength of the interaction. 
This is the ratio of a typical Coulomb potential to the thermal energy kT ||, 

(12) 



AnakT 



(with e 2 = Ana ~ 0.0917). In general < S > is a function of the density and temperature 
separately. However, if one ignores the relatively small effect of the screening length A e in 
Eq. @ then < S > only depends on T (and E). We have performed simulations for a pure 
56 Fe plasma at kT = 1 MeV. We assume results can be extrapolated to other compositions 
and temperatures by calculating the appropriate T. 

A least squares fit of our Monte Carlo results valid for all E v and 1 < Y < 150 is carried 
out. This fit is based on simulations for twelve values of T between 0.87 and 151.8. For a 
temperature of one MeV this corresponds to 56 Fe densities from 2 x 10 7 to 9 x 10 13 g/cm 3 . 
We approximate < S >, 

< S(E,T) >= l/[l + fflq>(-5>(r)#)], (13) 

i=0 

for 

E < E*(T) =3 + 4/r 1/2 . (14) 

While for E > E* we assume, 

<S(E,V)>=1. (15) 
The coefficient functions A(r), for i=3,4,5 and 6 are expanded in a power series in 

r l/2 

A(r) = Ai + A 2 r 1/2 + p i3 r + A 4 r 3/2 . (16) 

The coefficients are collected in table |. Finite size effects contaminate the Monte Carlo 
results for small E. Therefore we use RPA results for /3q, 

fa = ln[0.300/(0.300 + 3r)], (17) 

/3i = and = 6.667. 

The error in the fit is typically less then 0.01. Although, for very large T, < S > oscillates 
around one at large E. This oscillation is not reproduced by our fit and can lead to an error 
as large as 0.05. However, this only occurs at very high densities and is expected to have 
negligible impact on the dynamics. Again, the fit is valid for all neutrino energies and 
1 < T < 150. For smaller T a good estimate is provided by simply setting T = 1. (Note, 
here < S > is only important at very small neutrino energies.) Likewise, for T > 150 a 
reasonable estimate is provided by setting T = 150 (as long as the system is in the liquid 
phase). A solid is expected to form for T rs 180 This may be relevant for models of 
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type la supernovae JT3| . The very interesting problem of "Bragg diffraction" of neutrinos in 



a radioactive crystal remains to be investigated. Neutrino wave lengths can be comparable 
to the lattice spacing. 

We use this fit for < S > to calculate the mean free path of a neutrino in a plasma of 



ions, neutrons and electrons. For example, Cooperstein and Wambach [TJj] modeled matter 
at 10 12 g/cm 3 as consisting of X n = 6 percent free neutrons and 94 percent nuclei of average 
charge Z « 37 and average mass A m 97 at a temperature of 1.5 MeV. This is appropriate 
for the collapse phase of a supernova. We use this composition in calculating the mean 
free path. For simplicity, the composition and temperature are assumed not to change with 
density and we ignore the strong interactions between ions and or neutrons. 

The transport mean free path A is assumed dominated by elastic scattering off of nuclei 



and neutrons 15 



15 km , 10 \Ie\\2. . ,. C 1 . r , , 



A = (— =— ) [(l-X n )—<S>R e + X n (c n v 2 + 5c: 2 )]-\ (18) 

Pu x £j v > A 

Here p 12 is the density in units of 10 12 g/cm 3 , the weak couplings of a neutron c™ are given 



in ref. [TH] and R e is an additional correction factor that describes electron screening. This 



is calculated in ref. and can be approximated, see below Eq. 

,ctZ.. 1 



Re 



1 + 



V 



C / 1 + 2.5£ 2 A 2 



v e 



(19) 



Each ion has an electron cloud around it. Electron neutrinos or anti- neutrinos couple to this 
with strength, c e v = 2sin 2 0vy + |, while muon neutrinos do not see the electron cloud c e v ~ 
and thus R e ~ 1. 

The mean free path A is shown in Fig. 1. To our knowledge, almost all present supernova 
simulations use Eq. ( |i~8D with < S >= R e = 1. This leads to a very small mean free path 
(which traps neutrinos for densities of about p±2 = 0.5 and above). However, including 
< S > leads to a dramatic increase in A and to a large change in its density dependence. 
The rapid decrease of < S > with density can lead to a A which actually increases with 
density. Over a range of densities A for E v = 10 MeV is greater then 10 km. This is much 
larger then the unscreened A (~ 0.4 km at pu = 5). Finally, electron screening causes A for 
a v e to be about 15 percent larger then for a v^. 

Screening effects are even more important for lower E v . For example at 5 MeV, A is 
greater then 45 km even at p\2 = 10. This is larger then the size of the dense system (« 30 
km) so a neutrino-sphere may not form at all (for this energy). 

Figure 1 shows that A is larger then the size of the system for E v less then or equal 
to about 7.5 MeV. For E v between 7.5 and about 10 MeV the relatively large A will allow 
neutrinos to diffuse out of the system (in about a msec or less). These are main results of 
this paper. 

However, at E v = 20 MeV (or above) screening is reduced and the overall 1/E 2 scale of A 
is smaller so that the mean free path is significantly shorter. Screening is not very sensitive 
to temperature (as long as there are no large changes in composition). Changing T leads to 
a change in T see Eq. (|T2"D. However < S > is not a strong function of T. Likewise < S > 
is not very sensitive to the average Z of the material. Changes in the average A change a in 
Eq. ( |i~0|) and the overall factor C 2 /A in Eq. fllSP . Thus A decreases with increasing A. 
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The static structure factor S q describes the total strength to scatter from the medium. 
At high density, this may be dominated by the excitation of ion plasma oscillations rather 
then (quasi-) elastic scattering. These plasma osc. have energies to ~ (Z 2 e 2 n/M i ) 1 ^ 2 (with 
ion mass Mj) and may lead to a much larger net neutrino energy loss per v — A collision ||17|| . 
This can impact neutrino thermalization and heating. 

Screening effects will be all but absent after the supernova shock wave dissociates the 
nuclei. Then A will be relatively short because of scattering from large numbers of nearly 
free neutrons and protons. Thus, the neutrino opacity is small (because of screening) before 
the shock wave and large afterwards. This change in opacity could have important effects on 
the dynamics. Perhaps the situation is not unlike the photon opacity of the universe being 
large before and small after recombination. 

We now speculate on some of the implications of screening on supernova simulations. We 
emphasize that final conclusions await neutrino transport calculations. First, more neutrinos 
may escape allowing additional electron capture (as escape reduces the neutrino chemical 
potential). Electron capture reduces the pressure and the energy of the shock. 

For example a 20 MeV (or higher) electron capture neutrino produced at p\2 = 5 can 
have its energy reduced to about 10 MeV in a time of order 1/3 msec by electron scattering. 
The neutrino can then diffuse out of the core in about a msec. Alternatively the neutrino's 
energy can be reduced to 7.5 MeV (in about a msec) and then directly escape. 

Second, increased diffusion should raise the neutrino luminosity during the early stages 
of a supernova. [Perhaps, this could be observable if many prompt u e s are detected from 
a nearby supernova [0.] This may enhance the neutrino transport of energy to the shock. 
Screening has almost no effect on the opacity of low density matter or of the dissociated 
material after shock passage. Thus screening should not interfere with the ability of material 
near the shock to absorb energy from neutrinos. (Note, screening introduces a new density 
dependence to the neutrino interactions.) 

We have calculated the effect of ion screening on neutrino-nucleus elastic scattering. Our 
Monte Carlo results for the angle average of the static structure factor have been fitted to 
an analytic formula. This may allow the inclusion of screening in simulations. We find that 
the mean free path of a 10 MeV (or lower) neutrino is greatly increased. This could have 
important effects on the early stages of a supernova. 

We thank the Institut fur Theoretische Kernphysik for their kind hospitality. This re- 
search was supported in part by the US Department of Energy under Grant No. DE-FG02- 
87ER-40365 and by the Deutsche Forschungsgemeinschaft. 
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TABLES 

TABLE I. Parameters (3{j from a least squares fit of the angle averaged static structure factor 
< S >, see text. 



Coeff. j = 1 2 3 4 

-7.362056 0.5371365 -0.1078845 4.189612E-3 

/3 4 j 3.4489581 -0.40251656 9.0877878E-2 -3.4353581E-3 

5j -0.74128645 0.11019855 -2.5359361E-2 9.0487744E-4 

6j 5.9573285E-2 -1.0186552E-2 2.2791369E-3 -7.4614597E-5 
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FIG. 1. Neutrino transport mean free path vs. density. The solid lines include both ion 
< S 1 and electron R e 7^ 1 screening and are appropriate for z/ e , v e while the dashed lines 
for neglect electron screening R e = 1. Finally the dotted lines (used in most present supernova 
simulations) neglect all screening < S >= R e = 1. Top to bottom, the curves are for neutrino 
energies of E v = 5, 7.5, 10, 15, and 20 MeV. 
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